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ABSTRACT 

We derive power density spectra in the optical and X-ray band in the tirnescalc range 
from several years down to ~ a day. We suggest that the optical power density 
spectrum consists of two separate components: long timescale variations and short 
timescale variations, with the dividing timescale around 100 days. The shape of the 
short timescale component is similar to X-ray power density spectrum which is con- 
sistent with the interpretation of short timescale optical variations being caused by 
X-ray reprocessing. We show that the observed optical long timescale variability is 
consistent with thermal pulsations of the accretion disc. 

Key words: galaxies: active - galaxies:individual:NGC5548 - instabilities - X- 
rays: galaxies. 



1 INTRODUCTION 

Systematic multiwavelength study of radio quiet Active 
Galactic Nuclei (AGN) led to the emergence of the general 
picture of such a nucleus as a massive black hole surrounded 
by a relatively cool accretion disc of an outer radius below 
~ 1 pc and a hot optically thin plasma (see Mushotzky, 
Done & Pounds 1993). Comptonization of the soft photons 
by the plasma leads to emergence of strong X-ray emission 
from the nucleus. 

This compact unresolved region is surrounded by a 
dusty torus (see Antonucci 1993), marginally resolved in 
the IR at the distance scales below 100 pc in the case of 
the nearest AGN (see e.g. Greenhill et al. 1996 for water 
maser emission from NGC 1068) although not in the case of 
NGC 5548. 

The models of the compact innermost region have to 
reproduce the results of the detailed studies of the spectra 
of Seyfert 1 galaxies. X-ray observations were particularly 
successful in demonstrating the coexistence of the hot and 
cold gas close to the black hole through the detection of the 
reflection component (Pounds et al. 1990) and subsequent 
identification of the observed 6.4 keV iron line profile in a 
number of sources with the reflection from relativistic ac- 
cretion disc (Mushotzky et al. 1995 for IC 4329A and NGC 
5548, Tanaka et al. (1995) for MCG-6-30-15, Yaqoob et al. 
1995 for NGC 4151, Yaqoob et al. 1996 for NGC 7314, Nan- 
dra et al. 1997a for NGC 3516 and Weaver et al. 1997 for 
MCG -5-23-16; for a general discussion see Nandra et al. 
1997b). 

The shape of the high frequency tail of the X-ray spectra 



shows that the velocity distribution of electrons in the hot 
plasma is predominantly thermal (e.g. Gondek et al. 1996 
and Zdziarski et al. 1997 for averaged Seyfert 1 spectra, and 
Zdziarski, Johnson & Magdziarz 1996 for NGC 4151). 

However, the dynamics of this multi-phase medium 
feeding the black hole is still not understood although a 
number of attempts to construct a complete - usually sta- 
tionary - model were undertaken (e.g. Narayan, McClintock 
& Yi 1996, Witt, Czerny & Zycki 1997). The best way to 
learn which processes might dominate the flow is to study 
the variability of the nucleus since the variations can be 
directly related to dynamics whilst the spectra may be in- 
sensitive to dynamical parameters. 

Extensive studies of variability were done so far mostly 
for galactic X-ray sources, including black holes. In the case 
of AGN, most of the studies were constrained to the obser- 
vations of the spectral changes based on a few spectra and 
to the studies of the time delays between different spectral 
bands. Nevertheless those studies of a few monitored AGN 
brought important conclusions: the well correlated variabil- 
ity in the optical, UV, EUV and X-ray band and the lack of 
measurable delay between those frequency bands suggested 
that the optical/UV variability is caused by the reprocessing 
of the variable X-ray flux (e.g. Collin- Souffrin 1991, Gavel 
et al. 1992). 

However, detailed studies of that hypothesis were dif- 
ficult. Power spectra were constructed only for the X-ray 
data and early results based on long EXOSAT lightcurves 
showed that the PDS was a featureless power law without 
any characteristic timescale (McHardy & Czerny 1987 for 
NGC 5506, Lawrence et al. 1987 for NGC 4051). 
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The situation is changing during the last years. Due to 
the undertaken monitoring of several AGN more and more 
data in different wavebands are available, ft is therefore im- 
portant to study the physical constraints to the models im- 
posed by these data. 

ft is also important to compare the results of variability 
studies for AGN to what we know about the variability of 
galactic black holes. The question is the possibility of scaling 
of some, or most, properties of the flow with the mass of the 
black hole which may extend over orders of magnitudes, from 
10 10 of solar masses down to ten. If such a scaling exists for 
some properties, the difference in the timescales of the AGN 
and galactic black holes help to study them in a comple- 
mentary ways. Long evolution is easier to study for galactic 
sources, as the timescales of years for galactic sources corre- 
spond to millions of years for AGN. On the other hand the 
timescales shorter than the travel time through the horizon 
of a black hole, of order of tens or hundreds of seconds for 
AGN and of orders of tens of microseconds for a galactic 
black holes are easier to reach for AGN. 

In this paper we derive the power density spectrum of 
NGC 5548 in the optical band based on the data presented 
by Sergeev et al. (1997) including the data collected by AGN 
Watch team. We extend the PDS towards high frequen- 
cies using the nonlinear prediction method. We estimate the 
broad timescale PDS in X-ray band and we compare it to the 
optical one. We compare these distributions to the present 
knowledge of other AGN and of galactic X-ray sources. Fi- 
nally, we discuss the physical mechanisms of the observed 
optical variations. 



2 THE OBSERVATIONAL DATA FOR NGC 

5548 

The nearby Seyfert galaxy NGC 5548 (z=0.0174) is the 
most extensively observed AGN. The optical/UV variability 
of this object was followed systematically by International 
AGN Watch (Peterson et al. 1991; Peterson et al. 1992; Pe- 
terson et al. 1994; Korista et al. 1995) during the last five 
years. The X-ray variability was recently studied by Nan- 
dra & Pounds (1994) and Clavel et al. (1992) in the hard 
X-rays, by Done et al. (1995) in soft X-ray band, and new 
results are expected from XTE. The EUVE variability was 
detected by Marshall, Fruscione & Carone (1995) (see also 
Marshall et al. 1997). These data allowed to formulate first 
phenomenological models of the time behaviour of the flow 
in this object (e.g. Rokaki, Collin-Souffrin & Magnan 1993; 
Loska & Czerny 1997: Magdziarz et al. 1997). 

2.1 Optical data 

In this paper we use the data obtained by AGN Watch team 
(Peterson et al. 1994) together with data accumulated over 
the past ten years prior to the AGN Watch activity and 
made recently publicly available (Sergeev et al. 1997). 

The entire optical data set covers the period from 1977 
to 1993. Although the calibration of the data may not be 
very accurate (see Sergeev et al. 1997) the amplitude of the 
variability is large (close to factor 3) and small calibration 
errors of order of 3 percent cannot influence the general pic- 
ture. 



Basic properties of this lightcurve become apparent by 
simple visual inspection. The variations cover all timescales 
but the amplitude of variations is related to the timescale. 
In the shortest timescales, day-to-day changes do not exceed 
several percent. Within one year, large variations of order of 
40 percent peak-to-peak occurred. However, the lightcurve 
appears flat over the whole interval of 17 years, indicating 
possible deficit of fluctuations with longest time scales ex- 
ceeding several years. 

2.2 X-ray data 

The X-ray power density spectrum was already obtained 
by Papadakis & Lawrence (1993). It was derived only on 
the basis of relatively short observations made by EXOSAT 
satellite which covered only timescales below 0.3 d (more 
precisely, between 3 x 10 4 and 5 x 10 2 s), beyond the range 
covered by the optical data. The PDS was best described as 
a power law with additional QPO type feature at ~ 10 3 s. 
The slope of the aperiodic variability was not determined 
accurately even for the longest data set from July 1984 
(q = 1.99 ± 0.9). The variations in the 2-10 keV flux in 
this EXOSAT data of the length of 1.006 d reach the value 
0.42 x 10~ n erg/s/cm 2 (Kaastra & Barr 1989). 

We will try to estimate the X-ray PDS at longer 
timescales using the following information about the char- 
acteristic variability amplitudes. 

The variations observed by Ginga during the 42 days 
of monitoring reached the amplitude 1.58 x 10 _11 erg/s/cm 2 
(Clavel et al. 1992). The strong brightening of the source in 
the optical band in 1984 was accompanied by only a mod- 
erate brightening in X-rays showing that the amplitude of 
variations in timescales of ~ 6 years increases only up to 
1.73 x 10" n crg/s/cm 2 . 



3 RESULTS 

3.1 The optical power density spectrum for NGC 
5548 

We study the optical lightcurve of NGC 5548 using the 
Fourier analysis method. Although a number of other tech- 
niques were applied to study variability, determination of 
the power density spectrum still remains the basic approach 
even in the study of aperiodic variability (e.g. van der Klis 
1995). 

In the case of AGN only X-ray (and not the optical) 
variability was studied in such way, with exception of a BL 
Lac object OJ 287 (Silanpaa, et al. 1996), since the time 
extension and the coverage by optical data is usually not 
good enough to perform that kind of analysis. 

We calculate the power density spectrum (PDS) us- 
ing an algorithm of Deeming (1975) and applying sufficient 
sampling interval of 0.00001 dT . In order to reduce the 
noise and to enhance the signal we smooth the power spec- 
trum with a median filter and then bin it on the log scale. 
Following Korista et al. (1995) we adopt such normaliza- 
tion, that PDS is computed in units of flux squared times 
day, i.e. 10~ 30 (erg / cm 2 sAfd. At redshift z = 0.0174 the 
adopted unit of PDS corresponds to power in B band of 
30 x 10 18 L| Q d (c.f. Korista et al., 1995). No correction for 
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Figure 1. The smoothed and binned power density spectrum 
of NGC 5548 in the optical band plotted in the decimal log- 
log scales. On the horizontal axis frequency is in units of cy- 
cles per day. On the vertical axis power density is in units of 
10~ (erg /cm s A) d. Continuous line is a broken power law fit 
to the power spectrum, with the slope 1.17. 
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Figure 2. The smoothed and binned power density spectrum of 
NGC 5548 given as in Fig. 1, together with a direct fit (contin- 
uous line). Dashed lines show examples of numerical simulations 
probing the influence of the window function onto the fit. 

power leaks to high frequency induced by window function 
was applied. Simulations demonstrate that leaks are criti- 
cal for spectral indices approaching -2 (Green, McHardy & 
Lehto 1993). Major source of uncertainty in our power spec- 
trum are low frequency fluctuations for which statistics in 
our data is poor and which contribute large fraction of total 
power in the power spectrum. 

The resulting PDS is shown in Fig. 1. It is plotted 
in log-log axis usually adopted in X-ray variability analy- 
sis. Timescales shorter than ~ 30 days are dominated by 
white noise and better time resolution is required to study 
that range. We subtracted this constant component from 
the whole spectrum before plotting. The overall spectrum is 



dominated by clearly aperiodic signal. It can be described as 
a power law with index ~ — 1.17±0.04 flattening off at time 
scales longer than about 3 years. We show the fit of such a 
power density shape to the observed one (Fig.l, continuous 
line) . 

The above error estimates base on the internal scat- 
tering of the residuals from the fit. The influence of the 
window function itself on the shape of the spectrum was 
analysed by Papadakis & Lawrence (1995). They conclude 
that the power leak to high frequencies is negligible for in- 
clinations flatter than -2. Our own extensive simulations 
(Schwarzenberg-Czerny et al, in preparation) show that the 
actual external error of the slope may be as large as 0.2. 
However, the simulations confirm the existence of the break 
- its frequency is reproduced with accuracy logv = 3.2±0.3. 
Exemplary results of such simulations are shown in Fig. 2. 

At lowest frequencies our PDS is affected by smoothing. 
Hence we performed independent analysis of low-frequency 
variations. For this purpose we splitted data into two parts, 
namely AGN Watch data and pre-Watch data and then 
computed moments for each part separately. Variances in 
two parts are identical, corresponding to standard devia- 
tion of individual observation of a = 1.63. Corresponding 
averages, namely 8.61 and 10.05 differ by less than a stan- 
dard deviation. This may mean that PDS at lowest frequen- 
cies does not exceed average PDS value, in agreement with 
Fig.l. Note that because of correlation of individual observa- 
tions, we expect difference of ave rages co nsiderably in excess 
of white noise expectation a j yj (n — 1). This argument, as 
model independent, is stronger than any results from the 
simulations. Our simulations, however, support this conclu- 
sion (Schwarzenberg-Czerny et al, in preparation). 

Therefore, the flattening of the optical PDS at low fre- 
quencies is a real property of the data. However, we cannot 
tell from this analysis whether actually the change of slope 
is rapid, or gradual. The results of the nonlinear prediction 
method seem to support the second possibility (see Section 
3.2). 

Formally, our PDS indicates the presence of a period- 
icity, or quasi-periodicity, at a knee of the basic shape (i.e. 
~ 2.5 years), and its harmonics of 1.2 years. These features 
are still better seen before smoothing and binning of the 
PDS. However, the significance of that feature is doubtful 
for several reasons. On the one hand, the whole data cover 
only 7 periods for such a timescale. The power at this time 
scale comes from the two minima observed in 1990 and 1992 
as the PDS computed for data excluding this minima does 
not display this feature. Any process generating large fluc- 
tuations would produce such PDS features if by chance ex- 
actly two of them occurred during the observed interval. 
Let us assume for the moment that the probability of ob- 
serving minima in the lightcurve of NGC5548 is uniform in 
time. Then, the probability of finding two of them within 
2.5 years during the 17 years of observations is as high as 
2.5/17 - 0.15. 

3.2 The results of the nonlinear prediction 
method for NGC 5548 

The high frequency behaviour of the data is better analysed 
with the use of the nonlinear prediction method which is 
sensitive to the timescales as short as the typical separa- 
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Figure 3. The quality of prediction as a function of the prediction 
time-step for different timestep of the interpolated data: open 
squares (2.5 d, AGN Watch data), filled squares (5 d, AGN Watch 
data), open circles (10 d, AGN Watch), filled circles (10 d, whole 
data), stars (20 d, whole data), filled triangles (40 d, whole data) 
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Figure 4. The quality of prediction as a function of the predic- 
tion time-step: the combined picture. Data for prediction time 
below 25 d are from interpolation 2.5 d set, between 20 d and 50 
d from interpolation 5 d set, from 50 d to 100 d from interpola- 
tion 10 d set (for AGN Watch data and, separately, from whole 
data), between 100 d and 200 d from interpolation 20 d set and 
finally between 200 and 400 d from interpolation 40 d set. The 
correlations for the longest prediction time steps are close to zero. 
The dashed line shows the linear fit to the high frequency part 
well approximated by a power law. The slope of this power law 
equal to 0.74 corresponds to the PDS slope 1.74. 



tion between the measurements. This method, introduced 
by Tsonis & Eisner (1992) and successfully applied to X- 
ray data for AGN by Czerny & Lehto (1997) is based on the 
quantitative analysis of the decay of the correlation between 
any two short sequences in the data with an increase of the 
time separation of those sequences. It does not require any 
subtraction of the white noise. 

The correlation coefficient R is usually a decreasing 
function of the time separation T. If the signal is well rep- 
resented by a simple shot noise, its PDS is well described 
by a power law with an index a and at the same time the 
functional form of the dependence R(T) takes the shape 



1 - R(T) oc T a - 



(1) 



Therefore, in the case of non-periodic signal with ba- 
sically a power law shape of PDS the nonlinear prediction 
method can be used to determine independently the slope 
of the PDS. 

The method was developed for the equally spaced data 
with gaps so the original data had to be interpolated. 

To study the shortest timescales we used the data col- 
lected by AGN Watch team. The mean separation of the 
data points is 3.4 days, with sequencies sampled as densely 
as 1 day during the fifth year so we adopted the timestep 
2.5 d for interpolation. We included only these data points 
which were interpolated between the original data points 
separated not more than 1.6 times the adopted timestep. 

The resulting correlation coefficient R between the ac- 
tual and predicted values of the flux as a function of the 
prediction time T is shown in Fig. 3 (open squares). The 
log(l-R) vs. log(T) plot is generally of a power low shape 
although a conclusion about the presence of some curvature 
cannot be rejected. The slope of this power law is equal 0.74 
which translates into the slope of the PDS equal 1.74 be- 
tween the frequencies 0.04 and 0.4 dT 1 . 

The signal is clarly of stochastic nature and not od de- 
terministic chaos since the linear fit to log(l-R) vs. log(T) 
is much better (the correlation coefficient equal 0.95) than 
the linear fit to log(l-R) vs. T diagram (the correlation co- 
efficient equal 0.83). 

We repeated the analysis assuming the interpolation 
timestep equal 5 and 10 d (filled squares and open circles) 
which provided results for longer timescales. The slopes do 
not differ more than by 0.1 from the previous case although 
for longer timescales the slope is determined with larger er- 
ror. Also the statistics of the pairs for prediction becomes 
a problem. Therefore we repeated the analysis for 10 d in- 
terpolation timestep using the whole data, i.e. including the 
entire period of 17 years of observations (Fig. 3, filled cir- 
cles). We also attempted to detect the long timescale trends 
using 20 d and 40 d interpolation timestep but these results 
are already very close to R=0 line (no correlation) . 

In all computations the dimension of the string d (i.e. 
the length of the data sequence used for prediction; see Cz- 
erny & Lehto 1997) was obtained from the criterion that 
the correlation for prediction over a single timestep is the 
largest. In all cases d was equal 1. 

As the best representation of the data we choose the 
results obtained from 2.5 d interpolation step, supplemented 
at longer timescales by the computations made with the 
use of consecutive longer interpolation timesteps (Fig. 4). 
However, it is difficult to make any definite statement about 
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Figure 5. The X-ray power density spectrum in 2 - 10 keV in 
units 10 -22 (erg/ cm 2 /s) 2 d. The dashed line shows the solution 
with adopted low frequency slope equal zero and continuous line 
shows the solution with the steepest possible slope in low fre- 
quency band (see Section 3.3). Points mark the frequencies at 
which the amplitude was measured. 

the flattening of the curve at timescales about 30 d since the 
correlations at those timescales become very weak. 

To see better the comparison between the results based 
on nonlinear prediction method and on classical Fourier 
analysis we replot both on a single plot using the propor- 
tionality: 

PDSoc(l-R(T))T (2) 

and adjusting the normalisation to the Fourier results (Fig. 
6). The overlapping region is large enough to see that the 
results of the nonlinear prediction method extend nicely the 
Fourier derived PDS to high frequencies without any sud- 
den change of the slope. However, a better data would be 
necessary in order to see whether a step like feature at ~ 50 
days is present in the optical data or rather the PDS is well 
represented by a single component with a significant cur- 
vature. Further work of AGN Watch team may resolve this 
problem. 

3.3 X-ray power density spectrum 

In view of lack of a complete X-ray luminosity curve cover- 
ing all relevant time scales we consider maximum range of 
variability at a given timescale as more robust measure of 
variance. Given an estimate of the variance at several time 
scales, we proceed to estimation of power density function 
P(y) using Parseval equation, 

poo 

Var(y ) = / P{y)dv (3) 

which reduces the problem to the inversion of an integral 
equation. We defer derivation of Eq. (|^) and its inversion to 
Appendix. 



At the shortest timescales (below 0.68 day) we adopt the 
PDS slope equal 1.99 after Papadakis and Lawrence (1993). 
We normalize it on the basis of equation (2) using the vari- 
ability amplitude for the 2-10 keV flux in the EXOSAT data 
given by Kaastra & Barr (1989). 

We assumed that PDS between v — \d~ x and l/6yr~ 
may be represented by a broken power law with a break 
at v . The parameters of the power law are determined by 
inversion of Eq. (^), or more specifically of its particular 
form (see Appendix, Eq. ||). 

In the first solution we assumed c*2 = 0. Our solution 
yields break frequency i/rj — 6.0 x 10 -2 d -1 and the slope 
ctj = 1.71 (see Fig. 5, dashed line). The position of the v 
depends very strongly on the difference between the ampli- 
tudes at v = l/42d _1 and l/Qyr -1 and the last value is 
actually given by a single measurement so the derived value 
should be treated as an estimate of the order of magnitude. 
We also constructed solutions with arbitrary but nonvanish- 
ing value of 02- The solutions existed only if Q2 was smaller 
than 0.41. We show this limiting case in Fig. 5 (continuous 
line) . 

We are unable to exclude negative values of 0.2, for our 
method does not work well for non monotonic PDS. 

The two limiting cases for the solutions of the con- 
straints imposed on X-ray PDS well illustrate the overall 
character of the X-ray variability. Clearly, the PDS can- 
not be a single power law in the entire studied timescale 
range. However, we cannot tell whether the turnoff is sharp 
or rather we deal with a broad band curvature. 



4 DISCUSSION 

4.1 Character of the optical variability 

The optical variability is aperiodic and broad band covering 
all the available timescales from years down to a day. Present 
data do not allow to distinguish any QPO phenomena su- 
perimposed on the broad component although the traces of 
some substructure in the PDS seem to show up. 

The overall variability is dominated by the timescales 
around 3 years. The PDS would show it clearly if reploted 
as power per frequency decade i.e. log(Power x Frequency) 
versus log(Frequency) diagram, in analogy to popular rep- 
resentation of the radiation flux as log(Flux x Frequency) 
vs. log(Frequency) . In such a new plot the fiat part of the 
PDS change into an steeply increasing power law peaking at 
a turn off point at 3 years (See Fig. 7). 

The subsequent part, between the timescales from ~ 3 
years down to ~ 100 days looks almost flat since the power 
law index 1.17 is translated into just 0.17 in the new diagram 
so the PDS is really broad band. This part ends up possibly 
with a second peak at timescales of about 30 d. 

We also cannot reject the hypothesis that the PDS 
shows a continuous curvature with a single broad flat part 
covering the timescale range between ~ 30d and ~ 3y. 

At the shortest timescales (below 30 days) the decrease 
of the variability amplitude is much faster, with PDS index 
equal 1.74. 

The data may therefore represent a single variability 
mechanism operating most efficiently between the timescales 
of 3 years and 30 days. 
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Figure 6. The combined information about the power density 
spectrum of NGC 5548: stars are the optical PDS from Fig. 1, 
points of various shape mark the optical PDS derived from the 
nonlinear prediction method (Fig. 4) for timescales shorter than 
100 d, continuous line mark the X-ray PDS from Fig. 5. The 
normalization of the last two components were adjusted to match 
the high frequency tail of the first component. 

However, a flattening of the PDS at ~ 200 days with 
subsequent turn off at ~ 30 days may rather suggest that the 
optical PDS consists of two separate components. The long 
timescale component has a flat PDS above 3 years and turns 
down with an index somewhat steeper than the overall 1.17. 
At the shortest timescales there is an onset of another com- 
ponent which again is flat above 100 days, peaks somewhere 
at about 30 days and then turns off rather steeply. 

4.2 Comparison with X-ray variability of NGC 
5548 

The dominant timescale in X-rays is somewhere around 17 
days, although may be a factor of a few longer if the ampli- 
tude at 6 year timescale is underestimated. This is clearly 
different from the timescale of 3 years dominating the opti- 
cal variability. However, it coincides interestingly with some 
substructure present in the optical PDS at the timescale of 
~ 100 days. 

Particularly, if the optical PDS consists of two compo- 
nents, as suggested in Section 4.1, we can roughly identify 
the X-ray PDS shape with the short timescale component. 
In both cases the PDS is flat above 100 days, the turnover 
is around ~ 30 days (30 in Fourier approach, 30 in nonlin- 
ear prediction and 17 in X-ray data) and the slope after the 
turn-off is steep (-1.71 for the timescales of days in X-rays, 
-1.74 in optical band). 

In order to show better the comparison of the optical 
and X-ray PDS we placed the X-ray PDS on the same plot 
as the optical one (see Figs. 6 and 7). Since in that case 
units are clearly different we adopted the normalization of 
the X-ray PDS in that case as arbitrary and we adjusted it 



Figure 7. The same as Fig. 6 but after multiplication by the 
frequency to show the power per logarithmic unit of the frequency; 
the peak on this diagram shows the frequencies at which most of 
the power resides. 

to the high frequency tail of the optical PDS. The deficit of 
the power in X-ray PDS with respect to the optical PDS at 
longer timescales are now clearly seen. 

Therefore the data seem to suggest the possibility that 
on the timescales longer than ~ 100 days the optical variabil- 
ity is unrelated to X-ray variations while on the timescales 
below ~ 100 days the shape of the optical PDS is dominated 
by another component of the shape similar to X-ray PDS. 

4.3 Comparison with X-ray variability of AGN 

X-ray PDS for other AGN are determined more accurately. 
For us the most interesting cases are two Seyfert galaxies 
studied in the broad frequency range. 

NGC 5506 is a bright Seyfert 2 galaxy. The long history 
of the X-ray measurements of that source together with a 3 
day continuous observation by EXOSAT allowed to deter- 
mine its PDS in the entire range from ~ 10 years down to 
400 s (McHardy 1989). The shape of the PDS is basically 
that of a power law with index ~ 1.5, clearly flattening above 
the timescales of order of ~ 20 days. Its high frequency part 
between 10~ 5 — 10~ 3 Hz was actually better represented as 
a broken power law changing index from 1.0 to 1.9 towards 
high frequencies. 

The broad band X-ray PDS spectrum of NGC 5548 de- 
rived by us is therefore quite similar to NGC 5506. Even the 
position of the Vfi a t is similar which may not be surprising 
since the luminosities of the two objects are comparable. 

NGC 4151 is a nearby moderately bright Seyfert 1.5 
galaxy frequently observed in X-rays. The broad band PDS 
for that galaxy also shows clearly the flat part at timescales 
~ 100 days and a power law with index ~ 2.1 between 10 3 
and 5 x 10 4 s (Abraham & McHardy 1989; see also Pa- 
padakis & McHardy 1995). It is, however, impossible to tell 
on the basis of that data whether the flattening progresses 
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continuously towards low frequencies or a well defined knee 
is present somewhere at the timescale between a day and 
several days. 

Few other AGN were also studied basing on the long 
EXOSAT lightcurves which allowed to determine accurately 
the slope of the PDS for timescales between fraction of an 
hour and a day. Mean index of that slope determined from 
a sample of a few objects was 1.7 ± 0.5 (Green et al. 1993) 
or 1.55 ± 0.09 (Papadakis & Lawrence 1993). Studies based 
on Ginga data (e.g. McHardy 1989) as well as complemen- 
tary studies of the old EXOSAT lightcurves made using the 
non-linear prediction method (Czerny & Lehto 1997) showed 
that power law continues down to the timescales of 200 s. 
Recent study of the ASCA data for MCG 6-15-30 with the 
use of excess pair fraction method (Yaqoob et al. 1997) indi- 
cate that the power law with an index ~ 1.5 extends actually 
down to 20 s timescale. However, no significant change in the 
slope of the power law was noticed. 

Certain degree of success in modelling X-ray lightcurves 
of AGN was achieved using Linear State Space Models em- 
ploying first order autoregressive process AR[1] for mod- 
elling of dynamics (Konig & Timmer, 1997). This model 
accommodates flattening of the spectrum both at low and 
high frequencies. 

4.4 Possible mechanisms of optical variability 

The optical variability is dominated by the broad band of 
timescales starting from 3 years and extending down to 
shorter timescales with only slowly weakening power. 

This timescale is too short for typical viscous evolution 
which requires thousands of years (see e.g. Siemiginowska, 
Czerny & Kostyunin 1996). Therefore we look for explana- 
tion of the observed behaviour amongst the dynamical or 
thermal timescales. 

4-4-1 Dynamical instabilities 

The basic dynamical timescale is given by the Keplerian 
orbital motion. The mass of the black hole in NGC 5548 is 
most probably within the range 6 x 10 7 and 1.4 x 1O 8 M (see 
e.g. Loska & Czerny 1997, Kuraszkiewicz, Loska & Czerny 
1997 and the references therein). Therefore, we can adopt 
the value 1O 8 M0 as representative. 

The dynamical timescale is equal 3 years when the dis- 
tance x (expressed in units of the Schwarzschild radius) is 
equal 

x = 1600M 2/3 . (4) 

This distance translates into a distance of 16 light days for 
1O 8 M black hole. 

If the disc does not extend beyond that radius and it is 
clumpy or dynamically pulsating below it then the observed 
range of timescales is reproduced. 

There are no easy arguments why the disc should be 
truncated at that radius. However, low ionization broad 
emission lines (LIL), if formed at the disc surface, (e.g. 
Collin- Souffrin 1987) show the typical delay of 18 days (Pe- 
terson et al. 1994), i.e. about the required light travel time. 
Also the characteristic timescale for the X-ray variability 
and the short timescale optical component peaks somewhere 



around this value which may correspond to the light travel 
time across the hot plasma region but those two timescales 
can be also related in a natural way without any reference 
to the dynamical timescale in the disc (see Section 4.6). 

4-4-2 Thermal instabilities 

The thermal timescale at any given radius depends on the 
assumptions about the properties of an accretion disc. If 
we assume that the disc structure is well represented by 
the classical model of Shakura & Sunyaev (1973) then the 
thermal timescale is longer than the dynamical timescale 
by a factor a -1 , where a is the viscosity parameter. The 
location of the radius where the thermal timescale is equal 
to 3 years depends on the adopted value of the viscosity 
parameter. 

The only value suggested for that particular object is 
0.03 derived on the basis of the overall optical/UV/X-ray 
spectra (Kuraszkiewicz et al. 1997) within the frame of the 
accreting corona model. Adopting that value for a scaling 
we obtain the radius 

x = 160a^ 3 3 M 8 2/3 (5) 

which is at the distance of 1.6 light day from the central 
black hole of the mass 1O 8 M0. 

The inner parts of AGN accretion discs, within a frame 
of the a times the total pressure viscosity model, are ther- 
mally unstable. We can compare the radius given by equa- 
tion 4 with the radius where there is the onset of the insta- 
bility. 

Simple estimate of the transition radius from the radi- 
ation pressure to gas pressure domination gives (Shakura & 
Sunyaev 1973) 

x ab = 46a 2/ 2 3 1 M 8 - 2 / 3 M 16 / 2 2 5 1 (6) 

where we adopted the mean accretion rate 0.025 Mq/jt (af- 
ter Kuraszkiewicz et al. 1997) as a scaling factor. 

This value only weakly depends on the viscosity and 
is by a factor ~ 3 smaller than expected observationally. 
However, more accurate computations are encouraging. 

Our accretion disc model was computed taking into ac- 
count the vertical structure of the disc (Rozariska et al. 
1998). Free- free opacity was replaced by the appropriate full 
opacities (Kramers means) in the radiation transfer equa- 
tion. 

The disc is thermally (and viscously) unstable when the 
ratio (5 of the gas pressure to the total pressure is smaller 
than 2/5 (Shakura & Sunyaev 1976). Since we now compute 
the vertical structure (so the gas to the total pressure ratio 
is also a function of the distance from the equatorial plane) 
we estimate the tendency of the disc to develop the thermal 
instability by computing, at any radius, the ratio of the disc 
mass where /3 is smaller than 2/5 to the total mass. 

In Fig. 8 we show the plot of the fraction of the disc 
mass which is thermally unstable at a given radius. It is 
calculated as a ratio of the surface density of the gas with f3 
smaller than 2/5 to the total surface density. 

We see that the instability starts to be present at ~ 
200Rschw and at 150Rs c h w the entire disc become thermally 
unstable so the thermal instability of an a disc is a viable 
explanation of the shape of the long timescale optical PDS. 
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Figure 8. The fraction of the disc which is thermally unsta- 
ble plotted as a function of radius. Model parameters: M = 
10 s M e ,M = O.O25M /yr and a = 0.033. 



4.5 Possible mechanisms of X-ray variability 

The interesting conclusions about the nature of X-ray vari- 
ability were derived so far only for galactic sources. Recent 
XTE observations of Cyg X-l(Cui et al. 1997) applied to 
the concept of coherence function (Vaughan & Novak 1997) 
suggested that most of the time the hot medium itself is not 
strongly variable. The observed variability may be caused 
instead by scattering of the soft photons emitted close to 
the black hole by an extended medium (Kazanas, Hua & 
Titarchuk 1997). In this picture the timescale of the flatten- 
ing of the PDS measures directly the size of the comptoniz- 
ing region. 

In the case of NGC 5548 it would mean that the hot 
medium extends up to ~ 30 light days. 

One of the possibilities is that the Compton heated 
corona forms above an accretion disc (Begelman, McKee & 
Shields 1983). The outer radius of such a corona is located 
at ~ WOOORschw (i-e. 100 light days). The corona, however, 
scatters photons efficiently only up to ~ 10% of that distance 
(e.g. Ostriker, McKee & Klein 1991) which clearly makes an 
interesting explanation of the presence of the hot medium 
within the distance required by the shape of the X-ray PDS 
although it would be necessary to invoke some extra heating 
in excess of the Inverse Compton temperature. 

The fact that the representative emission radius of the 
Hp is comparable to the disc fraction covered by a corona 
may not be accidental as the presence of the corona scatter- 
ing the photons emitted closer to the black hole and redi- 
recting a part of them towards the disc surface enhances 
the emission of lines (see e.g. Kurpiewski, Kuraszkiewicz & 
Czerny 1997). 



4.6 Comparison with galactic X-ray sources 

Rapid aperiodic variability occurs in all types of X-ray bina- 
ries (for a review see van der Klis 1995). Since the observed 
X-ray flux from these objects is typically a few orders of 
magnitude higher than the observed X-ray flux from AGN 
the shape of the PDS for galactic sources could have been 
studied in much more detail than for AGN. Therefore, in or- 
der to make an adequate comparison we have to concentrate 
only on the basic charactersitics of galactic sources. 

The broad band PDS component consisting of a flat 
part below the frequency Vfi a t and a power law with index 
higher than 1.0 above it is characteristic for several different 
types of objects in different states. 

Black holes clearly display that kind of variability in 
their low (hard) states, with Cyg X-l being the best exam- 
ple. The value of Vfi a t is somewhere around 0.1 Hz but it 
may vary by a factor of a few between separate observa- 
tions (Belloni & Hasinger 1990). The power law part index 
is ~ 1.1 between Vfiat and ~ 3 Hz, and it ~ 2.0 at high 
frequencies. The value of Vfi at is 10 7 times higher than the 
vfiat °f tne optical PDS for NGC 5548, as it would be ex- 
pected from a simple scaling with the mass of the central 
black hole. This fact, together with the slope of the PDS 
suggests that the X-ray PDS for Cyg X-l in low state corre- 
sponds to the optical (and not X-ray!) PDS for NGC 5548. 
Similar effect can be present in the case of NGC 4151: X-ray 
PDS break at 14 d (Lawrence & McHardy 1995) is an or- 
der of magnitude too law to provide a perfect analogy with 
Cyg X-l. Unfortunately, optical PDS is not available for this 
source. 

What is more, the frequency of the change in the PDS 
slope of Cyg X-l is again by a factor 10 7 higher than Vfiat of 
the X-ray PDS of NGC 5548 thus indicating again a kind of 
substructure. Such a view is confirmed by recent high quality 
XTE observations of that source during the transition from 
the hard state to the soft state (Cui et al. 1997). In this 
case the power spectrum is clearly two component, with one 
component having Vfiat at 3 Hz (their Figs. 2 and 3) followed 
by a steep power law and second component, an underlying 
power law, relatively flat and without any clear flattening 
within the studied range of timescales (above 0.03 Hz). 

If the explanation of the shape of the optical and X-ray 
power density spectrum of NGC 5548 is correct (see Section 
4.4 and 4.5) we would expect the simple scaling with mass 
of the central black hole in the case of Vfi a t from X-ray PDS 
related to the size of the hot medium. The timescales related 
to accretion disc also scales roughly with the mass of the 
central black hole but they are at the same time considerably 
influenced by the luminosity to the Eddington luminosity 
ratio. 

The question, however, remains why the PDS of galacic 
Black Hole are similar to optical instead of an X-ray PDS of 
NGC 5548. In the soft X-ray range the answer may lie in the 
contamination of the spectrum by the direct contribution of 
the disc to the X-ray emission due to much higher temper- 
atures of GBH discs due to the scaling as M~ 1//4 with the 
mass of the black hole. However, in the case of hard X-rays 
the explanation should be different. It may be related to 
the fact that viscous timescales are relatively much shorter 
in galactic objects than in AGN. Therefore, in GBH both 
timescales can mix since the viscous timescale is only tens 



© 0000 RAS, MNRAS 000, 000-000 



Power density spectrum of NGC 5548 and the nature of its variability 9 



or hundreds times longer than thermal timescales. In AGN 
the viscous timescales are up to hundreds of thousands times 
longer than thermal timescales. From that point of view, the 
observational determination of Vfi at in all states is crucial. 

The same basic component also dominates the PDS of 
the island sources - low luminosity accreting neutron stars 
with negligible magnetic field (see van der Klis 1995). The 
value of Vfiat is in that case ~ 0.3 Hz, again as it would be 
expected from the scaling of the optical PDS of NGC 5548. 
It indicates that this kind of variability is not necessarily 
characteristic only for acreting black holes. 

Finally, similar component is also typical for disk-fed 
X-ray pulsars like Cen X-l, with the value of Vfi at again 
~ 0.3 Hz. 

However, in their high (soft) states black holes are char- 
acterized by a single power law PDS shape over several 
decades, with an index quite close to 1, as in the case o 
the optical PDS of NGC 5548 at intermediate frequencies. 
If any flattening is present it has to take place at frequencies 
lower than 0.001 Hz, at least in some sources. It may be too 
long for thermal timescales within the thermally unstable 
part of the disc. 

The qualitative changes are characteristic not only for 
black holes. Also accreting neutron stars change their vari- 
ability (and spectral) appearance with an increase of the 
accretion rate. Island sources change into banana sources 
and the PDS of upper banana sources - a single broad band 
power law with index only slightly larger than 1 - strongly 
resamble the soft state of Cyg X-l (van der Klis 1995). 

At this stage our comparison of variability of NGC 5548 
with galactic sources indicates both some similarities as well 
as differences and cannot be conclusive. 



5 CONCLUSIONS 

Presented results support the view that there are two clearly 
distinct physical mechanisms responsible for variability of 
NGC 5548. 

In the short timescales (below ~ 100 days) the vari- 
ability is caused by the Comptonization of the soft photons 
(emitted by the innermost part of accretion disc, i.e. very 
close to the black hole) by an extended hot medium, per- 
haps a corona. This scattering results both in variable X-ray 
flux and in variable optical flux due to absorption of a frac- 
tion of the X-ray flux by the accretion disc. The size of the 
hot gas region required by the shape of the power density 
spectrum is well estimated by the theory of the Compton 
heated corona. 

In the long timescales (above ~ 100 days) the optical 
variability is unrelated to X-rays. It is caused by thermal 
pulsations of the accretion disc. The size of the region of the 
disc, required by the shape of the power density spectrum, 
is consistent with the region dominated by the radiation 
pressure and therefore thermally unstable within the frame 
of a standard disc model of Shakura & Sunyeav (1973). 

The scaling properties of the variability with mass are 
difficult to confirm on a basis of a simple comparison of 
NGC 5548 with the galactic sources. The difficulties lie in 
the relative closeness of the viscous and thermal timescale in 
galactic sources (factor 0.01 instead of 0.00001 in AGN) and 
the contamination of X-rays in galactic sources by a direct 



emission of an accretion disc due to the systematic difference 
in the disc temperatures between AGN and GBH by two 
orders of magnitude. More research should be done on even 
approximate extension of the X-ray PDS of other AGN onto 
the timescales of years and on studies of the shape of PDS 
of dfferent galactic sources on the timescales longer than 
100 s in order to confirm or reject the presented physical 
interpretations of the variability. 
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Without loss of generality we may assume here that mean 
value was subtracted from data so that < / >= 0. Then 
LHS of Eq. (^) reduces to variance estimate < Var(v) >. 
Modulus of Fourier transform in RHS |JT(j/)| 2 may be ex- 
pressed by PDS P{y), but normalization requires care. We 
use Deeming (1975) Eq. (61) to ensure proper normaliza- 
tion. We take into account that, for considered St and shape 
of PDS, integration up to l/8t is equivalent to integration 
to infinity. Finally we obtain Eq. (H): 



T < Var(u) >= 2T / P(y')dv 



(3) 



In this paper we consider PDS consisting of two power laws, 
joined at v Q : 



P(u) = 



Po {V/Vo)-^- V < V a 



(4) 



where P , v , a\ and Q2 are constants. Then integration in 
Eq. (0) may be carried analytically yielding: 



< Var(v) >-- 



2PgV I J/_\ 1_ 

1 — ai \u } 



l — o>2 



l — a 2 



V > V 



+ < Var{v ) > v < v 



Equation (JE]) may be effectively solved for P(v) as long as 
vP(y) is monotonically decreasing with v. X-ray data at 
hand seem to be consistent with that. We substitute ob- 
served X-ray variances into Eq. (^) and solve for parameters 
a, P and v . 

This paper has been processed by the authors using the 
Blackwell Scientific Publications DTgX style file. 



APPENDIX 

Link between variance and power density spectrum (PDS) 
may be established using Parseval equation: 

f-OO 

\f(t)\ 2 dt= / \F(u)\ 2 dv (1) 



Let us consider model signal f(t) binned in time intervals 
St differing from null in the observed time interval T and 
vanishing elsewhere. Signal f(t) constitutes an example of a 
band limited signal. Its Fourier transform decreases outside 
frequency band < 1/T, 1/St >. Taking this into account and 
assuming reality of /, so that \J-( — v)\ 2 = |J-"(!/)| 2 one ob- 
tains approximation of Eq. ([j]) in the following approximate 
form 

/ ' \f(t)\ 2 dt^2 f ' \T{v)fdv (2) 

J-T/2 Jl/T 
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